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We report a detailed theoretical analysis of novel quadrupolar interactions observed between 
islands, which are disk-like inclusions of extra layers, floating in thin, freely suspended smectic 
C liquid crystal films. Strong tangential anchoring at the island boundaries result in a strength 
+ 1 chiral defect in each island and a companion —1 defect in the film, these forming a topological 
dipole. While islands of the same handedness form linear chains with the topological dipoles pointing 
in the same direction, as reported in the literature, islands with different handedness form compact 
quadrupolar structures with the associated dipoles pointing in opposite directions. The interaction 
between such heterochiral island-defect pairs is complex, with the defects moving to minimize the 
director field distortion as the distance between the islands changes. The details of the inter-island 
potential and the trajectories of the —1 defects depend strongly on the elastic anisotropy of the liquid 
crystal, which can be modified in the experiments by varying the material chirality of the liquid 
crystal. A Landau model that describes the energetics of freely mobile defects is solved numerically 
to find equilibrium configurations for a wide range of parameters. 

PACS numbers: 61.30.Jf, 61.30.Dk, 61.30.-v 



I. INTRODUCTION 



The interactions between colloidal particles, fluid 
drops, and phase-separated domains in liquid crystals 
(LCs) are mediated by the distortion of the ordered fluid 
and depend on the LC elastic constants as well as on the 
anchoring conditions at the inclusion boundaries. The 
tendency of such inclusions to form aggregates depends 
intimately on the properties of the inclusions as well as 
those of the LC, providing exquisite control of the self- 
assembled colloidal structures. A key feature of these 
novel systems is the presence of topological defects, which 
are singularities of the LC director field, close to the 
boundaries of the inclusions. In the strong anchoring 
regime, the director field is distorted from its equilib- 
rium uniform state near each inclusion and, in general, 
topological constraints lead to a distortion equivalent to 
that of a virtual +1 defect inside the inclusion. In or- 
der to match the uniform director field far away, one or 
more defects of opposite sign are nucleated in the liq- 
uid crystal matrix, outside the inclusion. The nature of 
these defects determines the symmetry of the director 
field around inclusions and their long-range interactions, 
for example, whether they are dipolar or quadrupolar. 
The short-range interactions are more difficult to predict 
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as the defects may move when the inclusions approach 
each other, leading to complex two-body interactions. 

When small spherical particles are suspended in three- 
dimensional (3D) nematics [l| .several beautiful colloid- 
defect structures are observed @, [1[ , including the "satel- 
lite" (a single —1 defect next to an inclusion with radial 
anchoring) 0, H, 0] , "Saturn-ring" (a defect-line around 
an inclusion with radial anchoring) 0], and "boojum" 
(a pair of strength 1/2 defects at the poles of inclusions 
with tangential anchoring) configurations Q. The direc- 
tor field around the inclusion-defect pair typically has 
dipolar or quadrupolar symmetry Q1, leading to a vari- 
ety of self-organized superstructures [ljj, El QjJ ■ 

In two dimensions (2D), the first experimental obser- 
vation of inclusions of nematic droplets in chiral smec- 
tic C films reported radial anchoring of the director at 
the droplet boundaries and single — 1 defects nearby [l3| . 
Later observations of nematic and isotropic inclusions in 
achiral Sm C smectic films reported tangential anchoring 
QUEUE Q], and either one pair of —1/2 defects or 
a single —1 defect at the inclusion boundaries. More re- 
cently, a temperature-dependent tangential to radial an- 
choring transition at inclusion boundaries in chiral smec- 
tics with small spontaneous polarization was observed 
[3 Ell i resulting in positional rearrangements of the de- 
fects along the inclusion boundaries. 

In principle, 2D systems of inclusions are simpler to 
model numerically than 3D systems as the computations 
required to calculate inclusion interactions coupled to 



mobile defects are considerably less demanding. Accu- 
rate simulations of a variety of 2D systems have been re- 
ported over a wide range of physical parameters. Impor- 
tant analytical developments based on techniques such 
as series expansions 1201. electromagnetic analogies [2lj |. 
and complex analysis [221 ] have provided a wealth of de- 
tailed information that complements the numerical anal- 
ysis and have been used to test the accuracy of the 
numerical results. The method of images was used re- 
cently in this context to calculate the long-range inter- 
actions between inclusions under a variety of conditions 
[23l ] and inclusion-defect pair configurations, at short and 
long range, have been calculated numerically with very 
high precision 0, HH, [26[ . The models reveal that the 
short-range interactions between quadrupolar inclusions 
are complex, with the defects moving as the distance 
between the inclusions is changed. When the islands 
are close together, the defects are shared equally by the 
two inclusions at equilibrium [2^ |. reminiscent of cova- 
lent atomic bonding where the electronic density peaks 
in the region between the atoms. (A similar 'bonding' 
mechanism where a circular defect line is shared by two 
inclusions, creating some spectacular equilibrium config- 
urations, was also predicted numerically and observed in 
3D nematics 0, [29| ) . 

In this article, we study the interactions between cir- 
cular islands, which are disk-like inclusions of extra lay- 
ers, in freely suspended Sm C films [3(| ■ These films 
are structured in 2D fluid layers, with the liquid crys- 
tal molecules aligned along an average direction that is 
tilted from the layer normal. The projection of the molec- 
ular tilt on the plane of the layers defines a 2D vector, 
the so-called c-director. The aim of this work is to pro- 
vide a quantitative description of recent observations of 
the interactions between topologically chiral islands in 
both chiral and achiral freely suspended Sm C films, in 
which the existence of an important class of inclusion in- 
teractions not studied previously in either 2D or 3D, was 
revealed, namely interactions that depend on inclusion 
chirality [3l[. The experiments showed that islands of- 
ten form linear chains with the dipoles aligned along the 
chains, similar to earlier observations, but that heterochi- 
ral island-defect pairs form a novel quadrupolar structure 
with the topological dipoles pointing in opposite direc- 
tions. Indeed, both dipolar chains and quadrupolar as- 
semblies may coexist on the same achiral Sm C film, the 
local organization being determined by the topological 
chirality of the individual islands. 

This paper is organized as follows. In Section II we 
briefly summarize the experimental results, highlighting 
those that will be compared with the numerical calcu- 
lations. In Section III we present a Landau free energy 
model that describes self-consistently the disposition of 
freely mobile defects in the film, and solve it numerically 
for a wide range of parameters in Section IV, where a 
detailed comparison with the observations is carried out. 
Finally, we discuss the results and summarize our con- 
clusions in Section V. 
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FIG. 1: Chiral islands in a freely-suspended smectic film 
of achiral liquid crystal. DRLM images of a smectic C 
film of racemic MX8068 showing (a) two islands containing 
+1 defects with the same handedness and (b) two islands 
with +1 defects of opposite handedness. The defects in (a) 
form a dipolar chain, while those in (b) form a topological 
quadrupole. The inset in (a) shows a single island with a 
+1 defect inside and a matching —1 defect in the background 
film, the pair forming a dipole. The laser illumination (shown 
by the arrow) is obliquely incident in the x-z plane, allowing 
us to differentiate between left- and right-handed islands by 
comparing the brightness of the brushes in their upper and 
lower halves. The equilibrium c-director field (h) is shown 
around two islands with (c) the same handedness and (d) op- 
posite handedness, the head of the h indicating the end of the 
molecule closer to the viewer. 



II. EXPERIMENTAL RESULTS 

Figure 1 illustrates characteristic equilibrium config- 
urations of island-defect pairs (with TV ~ 20 layers) 
interacting on a background film (N = 2) of racemic 
MX8068 [321 . Strong orientational anchoring at the is- 
land boundaries induces a +1 defect inside each island 
and —1 defects in the background film, each island- 
defect pair forming a topological dipole (inset, Fig. 1(a)). 
Such dipoles attract over long distances through an elas- 
tic dipole-dipole interaction and, when the islands are 
sufficiently close, reach an equilibrium separation where 
the elastic distortion is minimized. Two types of island- 
defect pairs are observed: either the pairs form linear 
chains with the dipoles pointing in the same direction 
along the chain, as in Fig. 1(a), or the dipoles point in 
opposite directions and the island-defect pairs form a 
quadrupolar structure, with the —1 defects placed sym- 
metrically above and below the line joining the islands, as 
in Fig. 1(b). A sketch of the equilibrium c-director field 
when two islands with the same handedness are close 
to each other is shown in Fig. 1(c), while when the is- 
lands are heterochiral, the quadrupolar structure shown 
in Fig. 1(d) is favorable [311 ]. 

The equilibrium center-to-center separation D of is- 
lands of radius R in dipolar chains was found experimen- 
tally to be close to 2\/2R, the value predicted theoret- 
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FIG. 2: Textures of heterochiral islands interacting on a film 
of 25% chirally doped MX8068. (a) The quadrupolar struc- 
ture is in equilibrium when the islands almost touch, (b) The 
equilibrium separation between —1 defects increases as the 
islands are separated using optical tweezers, (c) When the 
separation is sufficiently large, the quadrupolar symmetry is 
broken, (d) When the islands are forced even further apart, 
the quadrupole evolves into two separate dipoles. 



ically and confirmed in simulations using a free energy 
with a single elastic constant [2l|, [25|. Islands in the 
quadrupolar structure also exhibit a well-defined equilib- 
rium separation. When the island separation is altered 
using optical tweezers (33|, the vertical separation S of 
the —1 defects also changes. At sufficiently large island 
separation, the quadrupolar symmetry is broken and the 
structure evolves into two separate dipoles (3l| . 

The interaction between islands on an JV = 2 achi- 
ral Sm C film was found to be weak, with large orien- 
tational fluctuations of the c-director leading to large 
fluctuations in the position of the —1 defect. On chiral 
MX8068 Sm C* films, in contrast, director fluctuations 
are smaller and the interaction between islands is large 
enough to be measured using optical tweezers. On such 
chiral films, all of the islands are topologically homochi- 
ral and they form only dipolar chains. We propose that 
the reason for this bias is that the spontaneous polariza- 
tion of the chiral Sm C* phase prefers to point toward 
the liquid crystal-air interface at the edges of the islands, 
inducing only left-handed +1 defects. In films of racemic 
MX8068 doped with a small amount of chiral material, 
left-handed defects predominate but some right-handed 
ones are also observed. At 25% chiral fraction, for exam- 
ple, about 1% of the observed defects are right-handed. 
At this chiral fraction, the quadrupolar structure is found 
to be in equilibrium when the islands nearly touch (see 
Figure 2(a)). 

Optical tweezers may be used to displace the islands 
from their equilibrium separations. As illustrated in 
Figs. 2(b)-(d), the equilibrium positions of the —1 de- 
fects change simultaneously as the islands are separated. 
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FIG. 3: Equilibrium vertical separation S of the — f defects 
as a function of the island center-to-center separation D in 
the quadrupolar configuration regime, for racemic and 25% 
chirally doped films of MX8068, compared with the results 
of numerical calculations for systems with elastic anisotropies 
ft = 0.2 and k — 2.4. The straight lines are linear fits of the 
experimental data. The dashed lines are fits to the numeri- 
cal data shown here to guide the eye. The value of ft = 2.4 
was chosen to reproduce the slope of the linear fit for the chi- 
ral mixture. Note, however, that the theoretical results for 
this k underestimate the value of S/R at contact (D/R = 2). 
Fluctuations are expected to increase the average defect sep- 
aration, S, and account, at least in part, for this discrepancy. 
The average slope of the experimental data for the racemate 
is high, indicating very small values of ft, a limit which is 
numerically inaccessible in our simulations. 



In Figure 3 we show the experimentally measured vertical 
separation S of the — 1 defects as a function of the island 
center-to-center separation D for racemic and 25% chiral 
fraction films, plotted together with theoretical results to 
be described in the next section. 



III. ELASTIC FREE ENERGY OF SMECTIC C 
LAYERS 

The molecular orientation in smectic C free stand- 
ing films may be described by the c-director, a two- 
dimensional unit vector field c = (cosf?, sinf?) in the di- 
rection of the in-plane projection of the average molecular 
tilt. The total elastic free energy J 7 of a film of thickness 
I may be written as the sum of two terms: 

T = l I fedxdy + e [ f w ds, (1) 
Jn Jon 

where f2 is the area of the liquid crystal film, and dQ its 
boundary, including any inclusion boundaries that may 
be present. The second term accounts for anchoring at 
the boundaries and may be omitted in the strong anchor- 
ing regime. The first term describes the macroscopic free 
energy cost of elastic distortions of the smectic C director 
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field and is given by 

/ c = i(/C s (V-c) 2 +/C b (Vxc) 2 ) (2) 

where K, s and /Cj> are respectively the 2D splay and bend 
elastic constants. 

This macroscopic elastic free energy neglects variations 
of the degree of liquid crystalline order and it is inade- 
quate when disordering effects are important. This is 
the case when topological defects are present, as the or- 
der vanishes at the center of the defects. A better ap- 
proach, that is also more tractable numerically, is to use 
a mesoscopic description that takes into account varia- 
tions of the degree of order on the scale of the relevant 
correlation length. This is accomplished by theories of 
the Landau class that consider the free energy as a func- 
tional of a generalized order parameter, a quantity with 
the symmetry of the distortion field (XY or 2D vector 
in this case) that is non-zero in the ordered phase and 
vanishes in the disordered phase. This order parameter 
varies throughout a non-uniform system and is capable 
of describing the inhomogencitics in orientation and de- 
gree of order that occur at topological defects. If we 
replace the two-dimensional director field c = (c x ,c y ) by 
the complex order parameter <j> = <j) x +i<f>y = \(j>\(c x + ic v ) , 
we may write the bulk free energy (in reduced units [341 ]) 
as 

T = / f^dxdy , (3) 
Jq 

with energy density 

U = |0| 2 (|0| 2 -i) 

+ e [{d X <P X + dy(t>y) 2 + ~ dy^ X ) 2 } , (4) 

where £ is the correlation length (the characteristic size 
of the inhomogeneities), and k = ICt,/lC s the ratio of the 
elastic constants. This free energy includes all scalar in- 
variants of both the order parameter (up to fourth or- 
der in |0|) and its derivatives. The negative sign of the 
quadratic term in (f> indicates that the stable, uniform 
phase is ordered and in the absence of constraints the 
free energy is minimized by a uniform field of magnitude 
\4>\ = 1/V2. When the direction of 4> varies in space, its 
magnitude may also change in order to minimize the total 
free energy. At the defects, the order parameter </> varies 
rapidly, its magnitude vanishing within a distance of the 
order of the temperature-dependent correlation length £. 

In chiral films, the director field is additionally coupled 
to the free energy associated with the local space charge 
— V • P, where P is the spontaneous polarization, and is 
affected by screening from impurity ions dissolved in the 
film. We assume here that the effect of material chirality 
of the LC may be described by increasing the effective 
bend elastic constant, /Cf,, or, equivalently, by increasing 
the ratio of the elastic constants, k [35L l36j. 



IV. SIMULATION RESULTS 

In this section we present detailed numerical results 
describing the interaction between homochiral and hcte- 
rochiral islands for a range of parameters. In particular, 
we calculate the separation of the two —1 defects asso- 
ciated with a pair of heterochiral islands as a function 
of the island separation and compare this with the ex- 
perimental results to obtain an estimate of the effective 
elastic anisotropy of the experimental systems. 

We start by considering two islands with strong (fixed) 
tangential anchoring conditions, and the same handed- 
ness. For simplicity, we have taken \<p\ = l/\/2, the 
bulk equilibrium value of the order parameter for uniform 
alignment, at the island boundaries. At the borders of 
the domain, representing large distances, we have forced 
the director to be uniformly aligned along the y-direction. 
The correlation length that sets the scale for the size of 
the defects was taken to be £ = O.IR, a tenth the ra- 
dius of the circular islands, R, which were assumed to be 
monodisperse. This is admittedly an overestimate of the 
actual defect size but the major difficulty in the numeri- 
cal analysis stems from the disparity in the length scales 
R and £, and this choice is a compromise that speeds up 
the calculations and allows a systematic inve stig ation of a 
wide range of systems. As in previous work [251 ] . we have 
used finite elements with adaptive meshing to minimize 
the free energy. A first triangulation respecting the pre- 
defined boundaries is constructed. The order parameter 
<f> is specified at the vertices of the mesh and linearly in- 
terpolated within each triangle. The free energy is then 
minimized using standard methods, the Hessian of the 
solution at each iteration being used to generate a new 
mesh. If the solution presents stronger spatial variations 
(usually near the defects), a mesh refinement is required. 
On the contrary, if the solution is almost constant over a 
large region (far from the islands, for example), a coars- 
ening of the mesh is possible. In typical calculations, 
convergence is obtained after only two mesh adaptations, 
corresponding to final meshes with 6 x 10 3 points, span- 
ning a region of 15R x 15i?, and minimal mesh sizes of 
10 _4 i? close to the defects. The free energy is obtained 
with a relative accuracy of 1CP 4 . 

We have calculated the free energy of a pair of ho- 
mochiral islands forming a chain with their satellite 
— 1 defects as a function of the island center-to-center 
separation D. In order to investigate the effects of the 
elastic anisotropy, k was varied over a wide range. We see 
from Fig. 4(a) that the interaction between two islands is 
attractive at long distances, exhibits a well defined mini- 
mum, and is repulsive at small distances. For k = 1, the 
equilibrium separation D m i n = 2.80R is slightly different 
from that obtained in our previous work, D m i n = 2.82i?, 
based on the Frank elastic free energy ([T][with no order 
parameter variation and pinned defects [25| • This differ- 
ence is a consequence of the larger value used for £ here. 
Nevertheless, in order to carry out extensive calculations 
and explore a reasonable range of parameters, we decided 
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FIG. 4: (a) Total free energy as a function of homochiral is- 
land separation D, for different k. When k = Kf,/K 3 = 1, the 
free energy exhibits a well-defined minimum at a separation 
2y/2R. (b) Equilib rium separation D m m as a function 
of the elastic anisotropy n of the system. 
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FIG. 5: Equilibrium director configurations for different het- 
erochiral island separations with k = 1. (a) D = 2.2R: the 
defects are equally shared by the islands, in a symmetric, 
quadrupolar configuration, (b) D = 3.07?: the symmetry is 
broken; here the defects are found in a metastable configura- 
tion where they are both closer to one island, (c) D = 3.57?: 
each defect follows a different island, (d) D = 4.07?: even at 
this separation, the positions of the defects are still strongly 
influenced by the presence of the other island-defect pair, 
with configurations typical of isolated dipoles occurring at 
larger separations. 



to keep £ = 0.1R, as the discrepancy in theoretical equi- 
librium distances obtained is small compared with the 
uncertainty in the experimental data. The equilibrium 
separation D m i n as a function of k is plotted in Fig. 4(b), 
showing the separation increasing sharply as the bend 
elasticity is reduced below (n — > 0) but approximately 
constant for «; > 1. 

We next considered two islands with opposite handed- 
ness, with the defects located between them as observed 
in the experiments and shown in Fig. [2] In this case, 
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FIG. 6: Equilibrium director configurations for different hete- 
rochiral island separations with n — 0.2. The order parameter 
around the —1 defects has no longer a circular symmetry as 
in Fig. [S] This is reflected in the narrowness of the white 
brushes that oome out of the —1 defects, (a) D = 2.2R. (b) 
D = 3.07?. (c) D = 3.57?. (d) D = 4.07?. 



the defects strongly influence the interaction between the 
islands, and they move considerably as the distance be- 
tween the islands changes. As a result, the number of 
mesh adaptations had to be increased to three or four 
in order to obtain reasonably accurate results (indicated 
by the smoothness of the free energy plots, for exam- 
ple). Figures O and [5] illustrate model director configura- 
tions for different island separations D for systems char- 
acterized by two different elastic anisotropies, k = 1 and 
k = 0.2. When the islands are very close together, the 
defects are located symmetrically in a quadrupolar struc- 
ture, being shared equally by the two islands (Figs.[5^a), 
|U(a,b)). Above a certain island separation, this symme- 
try is broken and each defect becomes associated pri- 
marily with just one of the islands. We also sometimes 
observe configurations in which both defects are closer 
to the same island (see Fig. EJb)), although these are 
metastable. The defects influence each other over large 
distances (Figs. EJc,d), EHc,d)), their relative positions 
with respect to their parent islands being very different 
from an isolated island-defect configuration. When the 
elastic anisotropy of the liquid crystal is large, the defect 
brushes become sharper ,371] ■ as shown in Fig. [H]for the 
case JC S = 5/Cf, (n = 0.2). The interaction between the 
islands is, however, qualitatively similar to the isotropic 
elasticity case. 

The effect of elastic anisotropy on the interaction be- 
tween two heterochiral islands was investigated by calcu- 
lating the free energy as a function of the ccntcr-to-ccnter 
island separation, D, for different values of k. The results 
are plotted in Fig. E£a). For small k, the potential is sim- 
ilar to that observed with homochiral islands (Fig. SJa)) 
but the equilibrium separation is smaller and decreases 
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FIG. 7: (a) Total free energy as a function of heterochiral 
island separation D, for k = ICb/IC s = 1, 2, and 3. All systems 
have a well defined equilibrium separation Z? m i n that is smaller 
than in the homochiral case and decreases monotonically with 
increasing k (b). 

monotonically with shown in Fig. [7(b) . The islands 
are much closer together than in the homochiral case and 
when k is large, the islands almost touch each other, as 
observed in the experiments [3l| . 

We next studied the geometric evolution of heterochi- 
ral island-defect pairs as the island separation was var- 
ied. In particular, we tracked the vertical separation of 
the — 1 defects as a function of the distance between is- 
lands in order to compare it with the experimental results 
and to estimate the effective elastic anisotropy of the ex- 
perimental systems. 

Determining the trajectories of the —1 defects as a 
function of island separation is computationally very de- 
manding, as an accurate determination of the position of 
the defects requires a large number of mesh adaptations, 
and is particularly difficult at the quadrupolar symme- 
try breaking transition, an effect we now describe. At 
small separations, the islands share the defects in a sym- 
metric, quadrupolar configuration. In this regime, the 
defects are along the vertical bisector of the islands, so 
their horizontal coordinate = (see Fig. [5(a)) but 
at a critical island separation, D c , the quadrupolar sym- 
metry is broken and Xd becomes non-zero. One of the 
defects moves to the right, acquiring a positive Xd offset, 
while the other moves symmetrically to the left and has 
negative Xd- This behavior is qualitatively the same for 
all values of the elastic anisotropy k that were investi- 
gated, but the island separation at which the transition 
occurs depends on n. Fig. [5(b) shows the critical island 
separation D c as a function of k. The trend observed in 
D c confirms that the quadrupolar symmetry is broken in 
a continuous transition, or bifurcation, that occurs at a 
precise value of the island separation. 

Similarly, the vertical displacement of the defects yd 
was determined as a function of D, with results shown 
in Fig. [9(a) for several values of k. For moderately 
anisotropic systems (k = 1.4 and 2.4), the vertical dis- 
placement of the defects first increases, approximately 
linearly, as the island separation increases. When Xd 
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FIG. 8: (a) Computed horizontal position of the right defect 
X d £LS £L function of island separation D, for k = 1.8, 3.0, 
and 4.0. The reference system has its origin at the mid-point 
between the islands, (b) The quadrupolar symmetry of the 
island-defect pair is generally broken at a critical island sepa- 
ration D c , shown here as a function of k. The noisiness of the 
data when k < 1 reflects the numerical difficulties encoun- 
tered in modeling this elasticity regime. 



bifurcates at D c w 2.5R, however, the behavior of yd 
changes, increasing to a maximum and then decreasing 
slowly, until it reaches the isolated island-defect value, 
yd = 0, at large separation. These results illustrate how 
strongly the island-defect pairs interact even at quite 
large distances: for D ~ 5R, the vertical displacement 
of the defects is still very large, yd ~ R- The initial slope 
of the yd vs. D curve, which is a measure of the suscep- 
tibility of yd to fluctuations in the island separation, is 
plotted in Fig. 9(b) as a function of k. For values smaller 
than k = 0.2, the numerical results become unreliable: 
although the energy can still be obtained accurately in 
this regime, the equilibrium defect positions are sensi- 
tive to numerical noise and can not be uniquely defined. 
As the islands separate, the position of the defects may 
vary substantially, and due to the high degree of accuracy 
required for small k, our method is not able to resolve 
them. 

The numerical vertical separation of the defects, S = 
2yd, were compared with the experimental results in 
Fig. [3] The value of n — 2.4 was chosen to reproduce 
the slope of the linear fit for the chiral mixture. Note, 
however, that the theoretical results for this k underes- 
timate the value of S/R at contact (D/R = 2). Keeping 
the same k, we decreased the correlation length £ by a 
factor of ten. However, we still obtained a smaller value 
of S/R at contact. Fluctuations are expected to increase 
the average defect separation, S, and account, at least 
in part, for this discrepancy [23j|. The average slope of 
the experimental data for the racemate is high, indicat- 
ing very small values of k, a limit which is numerically 
inaccessible in our simulations. 

In principle, the susceptibility of the vertical displace- 
ment of the —1 defects to perturbations of the islands 
from their equilibrium separation plotted in Fig. [H(a) 
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FIG. 9: (a) Computed vertical displacement yd of the upper 
— 1 defect as a function of island separation D, for re = 1.4, 2.4, 
and 4.5 [j/d is half the vertical defect separation S plotted in 
Fig- El] For small D (close to equilibrium), the displacement 
of the defect increases linearly with a well-defined slope m, 
shown in (b) as a function of the elastic anisotropy re. 



should allow an estimate of re in the experiments. The 
calculations suggest a strong variation of this suscepti- 
bility as the elastic anisotropy re changes, as shown in 
Fig. EJb), but we have thus far been unable to reconcile 
the predicted behavior quantitatively with the experi- 
mental data. 



V. CONCLUSION 

We have described numerical simulations of the spon- 
taneous quadrupolar symmetry breaking of the director 
field around heterochiral island-defect pairs in Sm C liq- 
uid crystal films. In these films, the c-director is tan- 
gentially anchored at the island boundaries, inducing a 
+1 defect inside each island and a —1 defect outside in 
the background film, the island-defect pair forming a 
dipole. Such dipoles attract each other over long dis- 
tances through an elastic dipole-dipole interaction and, 
when the islands are sufficiently close, reach an equi- 
librium separation where the overall elastic distortion is 
minimized. 

Two types of island-defect pairs observed in the exper- 
iments were considered in the theoretical analysis. Either 
the islands have the same handedness, in which case they 
form linear chains with the dipoles pointing in the same 
direction along the chain, or the islands have different 
handedness and the dipoles point in opposite directions. 
In the latter case, at short distances, the island-defect 
pairs form a quadrupolar structure, with the —1 defects 
placed symmetrically above and below the line joining 
the islands. 

The simulations show that interactions between two 
heterochiral islands occur over a much longer range than 
homochiral island-defect pairs. Even when the islands 
are separated by up to five times their radius, the interac- 
tion between heterochiral island-defect pairs is still quite 



different from the long range dipole-dipole interaction ex- 
pected in the linear asymptotic regime, and an island- 
defect configuration significantly different from that ex- 
pected for isolated dipolar pairs. 

We have observed, both experimentally and numeri- 
cally, a transition where the quadrupolar symmetry of 
the heterochiral island system is broken. This occurs at 
a well-defined critical island separation that depends on 
the effective elastic anisotropy of the liquid crystal and 
the simulations suggest that the quadrupolar symmetry 
is broken continuously. For each system, characterized by 
an elastic anisotropy re, a well-defined critical island sepa- 
ration was found where the defects are displaced symmet- 
rically from the vertical mid-plane bisecting the islands. 
The island separation at which this transition occurs de- 
creases smoothly as re increases (at least for re > 0.5), 
giving further support to the existence of a bifurcation 
that delimits the quadrupolar configuration regime. The 
numerical calculations of the defect positions are compu- 
tationally demanding, leading to somewhat noisy results 
for systems with small re < 1. Furthermore, at the tran- 
sition large thermal fluctuations of the director field are 
expected and were observed experimentally in racemic 
systems. Including their effects in the model will affect 
the location (and maybe even the nature) of the bifurca- 
tion predicted by our mean- field analysis. 

At small island separations, the computed vertical dis- 
placement of the defects was found to increase linearly, in 
general agreement with experiment. The scatter in the 
experimental data strongly suggests that thermal fluctu- 
ations in the position of the defects are important even 
in the quadrupolar configuration regime and that these 
fluctuations should be taken into account in a full de- 
scription of the symmetry breaking transition. The sim- 
ulations suggest that the initial slope of the S vs. D curve 
may be used to estimate the effective elastic anisotropy 
re = Kb/K-g of experimental systems. We have found, 
however, that a good fit to the slope underestimates the 
defect separation at contact. It is likely that thermal 
fluctuations, which increase the average separation be- 
tween defects, are responsible for this discrepancy. Fur- 
thermore, the polydispersity of the islands may also af- 
fect the simulation results reported here. Heterochiral 
islands of different sizes do not exhibit perfect quadrupo- 
lar symmetry and this may result in a different behaviour 
of the experimental defect separation. Islands with dif- 
ferent sizes and the effect of thermal fluctuations will be 
investigated in future work. 
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